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Abstract 

We consider the one-loop factorization of the simplest twist-three process: inclusive deep-inelastic 
scattering of longitudinally-polarized leptons on a transversely-polarized nucleon target. By study- 
ing the Compton amplitudes for certain quark and gluon states at one loop, we find the coefficient 
functions for the non-singlet twist-three distributions in the factorization formula of g2( x B,Q 2 )- 
The result marks the first step towards a next-to-leading order (NLO) formalism for this transverse- 
CN ■ spin-dependent structure function of the nucleon. 
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Deep-inelastic scattering (DIS) of leptons on the nucleon is a time-honored example of 
the success of perturbative quantum chromodynamics (PQCD)]I[. The factorization for- 
mulae for the leading structure functions Fi(xb,Q 2 ) and gi(xs,Q 2 ), augmented by the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations for parton distri- 
butions 0, can describe the available DIS data collected over the last 30 years exceedingly 
well. Although the same formalism is believed to work for the so-called higher-twist struc- 
ture functions ||, e.g. g2(xs,Q 2 ) and Fl(xb,Q 2 ), which contribute to physical observables 
down by powers of the hard momentum Q, there are few detailed studies of them in the 
literature beyond the tree level. The QCD radiative corrections to g2(xs,Q 2 ) need be in- 
vestigated as accurate data have recently been taken j3J and more data will be available in 
the future ||. 

In this paper, we report a one-loop study of inclusive deep-inelastic scattering of 
longitudinally-polarized leptons (e.g. electrons) on a transversely-polarized nucleon tar- 
get ||. The subject was first investigated in the context of single parton scattering in Ref. 
0, and studies along the same line have continued in the literature ||. However, due to the 



subtlety of the twist-three process || [LC| [TlJ], those results are sensitive to the treatment of 
quark masses and are incomplete in the context of QCD fatorization. Indeed, even at tree 
level one must go beyond the single quark process to derive the correct ^2(^5, Q 2 ) expres- 



sion in terms of the parton distributions [12], [L3|, [L4|]. When loop corrections are included, 
one needs a general strategy to systematically calculate their contribution to higher twist 
processes. 

For a transversely polarized nucleon of four-momentum P^ and polarization vector S±, 
the hadron tensor = ± J e iq< (PS\[J tl (^), J u (0)}PS) can be expressed as 

W r = _ ie »v«p qaS±/3 l ( gi (x B ,Q 2 ) + g2(x B ,Q 2 )) , (1) 

where g M is the photon four-momentum, Q 2 = —q 2 , v = q- P, xb = Q 2 /{2v), and e 0123 = +1. 
J M is the electromagnetic current of the quarks. Thus it is the combination qt{xb,Q 2 ) = 
g±(xB, Q 2 )+g2(%B, Q 2 ) that naturally appears in the 1/Q-suppressed transverse polarization 
asymmetry. In this paper, we concentrate only on the non-singlet part of Qt{xb,Q 2 )', the 
singlet case is left to a separate publication. We seek a factorization formula for gr{xB,Q 2 ) to 
one-loop order in terms of the perturbative coefficient functions C{(x, y) and the generalized 
parton distributions (correlations) Ki(x,y) with two light-cone (or Feynman) variables x 
and y. The starting point is one-loop forward virtual-photon Compton scattering off a few 
"on-shell" quark and gluon states. ^From the scattering amplitudes, we examine the validity 
of infrared factorization and extract the one-loop coefficient functions Ci(x,y). 

We begin by outlining a general approach to the factorization of higher-twist observables, 
generalizing the method of Ref. ]TS] for the twist-two structure functions. For deep-inelastic 
scattering, it is convenient to consider first the forward Compton amplitude 

= ij d 4 £e iq< (PS\TJ^)J u (0)\PS) . (2) 

The spin- dependent part of the Compton amplitude (antisymmetric in /z and v) defines two 
invariant amplitudes Si^xb, Q 2 ), 

V = -ie^q a ( S ^^ S , + S J^1 {vSp - q ■ SPp)) • (3) 
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FIG. 1: Schematic diagram for Compton scattering in terms of quark and gluon scattering ampli- 
tudes and their correlation functions in the nucleon. 



According to the optical theorem, the imaginary part of the S\ y 2 amplitudes (divided by 
2it) is just the nucleon spin structure functions gi )2 . Hence, factorization of St{xb,Q 2 ) = 
Si(xb, Q 2 ) + S 2 (xb, Q 2 ) naturally leads to factorization of the structure function gr(xB, Q 2 )- 
The former has a straightforward Feynman-Dyson perturbative expansion, making it easily 
handled. 



According to Ref. [15), the nucleon Compton amplitude T^ u can be expressed as a sum 
of terms that are convolutions of quark-gluon Compton scattering amplitudes Mf v and the 
bare quark-gluon correlation functions I\s m the nucleon, as shown schematically in Fig. 1, 

TV = Y d M? ®YiB . (4) 

i 

Implicitly involved in the convolution are integration over the intermediate quark-gluon 
four- momenta kj and summation over the spin and color indices. M^ u , without external 
quark-gluon legs and self-energies, is the sum of a complete set of Feynman diagrams for 
Compton scattering. These diagrams are calculated in unrenormalized perturbation theory 
in the sense that all parameters in the expansion are bare. We use dimensional regularization 
to regularize both infrared and ultraviolet divergences in the diagrams. 

The pair of light-cone four-vectors p» = A(l, 0, 0, l)/y/2 and n» = (1, 0, 0, -l)/(>/2A) 
define a collinear basis in which the nucleon and photon momenta are written P M = p^ + 
M 2 n tl /2 and = —xsP^+vn", respectively, where A is an arbitrary dimensionful parameter. 
Intermediate quark and gluon momenta can also be expressed in this basis: 

= (kj ■ n)p» + {kj • p)n" + k f , (5) 

and the full result for M- IU can be expanded about kj = (kj ■ n)p^. The leading term in 
this expansion can be interpreted as scattering of collinear partons with Feynman momen- 
tum fractions Xj = kj ■ n. Because the partons are on-shell, M[ u can be viewed as the 
parton scattering S-matrix element. In principle, one must multiply by parton wave func- 
tion renormalization factors to get the proper S-matrix element. However, in dimensional 
regularization, the absence of a physical scale at the massless poles of the quark and gluon 
propagators reduces these contributions to unity. Subleading terms in the expansion of Mf 
are parton scattering S-matrix elements with insertions of certain vertices associated with 
powers of kj- and kj(= p ■ kj). For instance, with one power of quark momentum kj, the 
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subleading term is calculated with one insertion of the vector vertex to one of the quark 
propagators. Because S-matrix elements and their relatives are gauge invariant, one may 
choose any gauge for the internal gluon propagators in Mjf . We use Feynman's choice in 
our calculation. 

After the collinear expansion and integration over the quark-gluon four momenta, the 
convolution of Eq. (f|) involves integrating over the parton Feynman variables Xj = (kj ■ 
n). The correlation functions Y^b are matrix elements of gauge- invariant nonlocal quark 
and gluon operators in one-to-one correspondence with the external parton states in M^ v . 
In particular, when M^ u contains momentum-related vertex insertions, the quark and/or 
gluon fields in Fib appear with partial derivatives. The collinear expansion results in all 
QCD fields separated in spacetime along the light-cone direction n M . As in the leading- 
twist case, expanding the gluon polarization indices and summing over all contributions 
from longitudinally-polarized gluons generates straightline gauge-links which connect fields 
at separate spacetime points. To simplify the derivation of a factorization formula, we 
use the light-cone gauge (A • n — A + = 0), or link-free, expression for T iB . This gauge 
choice allows one to focus on the physical partons without worrying about the effects of the 
longitudinally polarized gluons. The gauge-invariant form of the factorization is recovered 
simply by imposing gauge invariance on the final form of parton correlations. 

As we have argued above, M^ v is ultraviolet finite because the on-shell wave function 
renormalization is trivial in dimensional regularization. Nevertheless, due to the massless 
on-shell external states, M% v has infrared divergences showing up as 1/e poles. These 
divergences may be factorized in the perturbative sense 

Mr = cr ® ^ , (6) 

where is the finite coefficient function and Pj contains only the 1/e poles. On the other 
hand, the infrared-finite quantities T^b contain ultraviolet divergences which also show up 
as 1/e poles. When the infrared poles in Pj cancel all the ultraviolet poles in T^ v is said 
to be factorizable. The product PjYiB defines the renormalized parton correlation functions 
Tj. The final factorization formula for the Compton amplitude is then 

= J2 c r ® T i > (7) 

i 

where is a well-defined perturbation series in a s and Tj is a finite nonperturbative 
distribution. 

We apply the above discussion to Compton scattering on a transversely polarized nucleon. 
As shown in Fig. 2, the non-singlet twist-three process involves two possible intermediate 
parton states in the light-cone gauge. First is the two-quark state with the transverse 
momentum flowing through , as shown in Fig. 2a. After Taylor-expansion, we keep only 
the contribution with exactly one insertion of the 27" vertex. The corresponding correlation 
function is 

r>(x )2/ ) = J ^e lXx e^y- x \PS\midl*P(\n)\PS) 

= 8(x - y) / —e ax (PS\m^(\n)\PS) , (8) 

where the Dirac and color indices on quark fields are open and a is perpendicular to p M and 
n M . The second intermediate parton state involves two quark lines and one gluon line, as 
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(a) (b) 

FIG. 2: Parton intermediate states appearing in the twist-three Compton amplitude. 

shown in Fig. 2b. If we use the Feynman rule it a ^f_ for the gluon attachment to M^ q) the 
corresponding correlation function is 

^ gqB (x,y) = J ^^e^^PSim^Al^i^PS) . (9) 

We remind the reader that all fields and the couplings here are bare. To decouple the spin 
and color indices in M^ v and IYb, we write explicitly 

r qgqB( X >y) = { S ±^ ^l(qgq)B (x, y) + iT£ ^2(qgq)B {x, y)) , ( 10 ) 

where N c is the number of colors, Cp = (N% — 1)/(2N C ), and T a = e a ^ lS Sx/3Pjns. In 
dimensional regularization, 75 must be defined explicitly. Here, we follow t' Hooft and 
Veltman's convention. 

With appropriate insertions of light-cone gauge links, which can be generated by sum- 
ming over intermediate states with additional longitudinally-polarized gluons, we define the 
following gauge-invariant parton correlations (TT2], [H]: 

K iB (x, y) = T i(qq)B (x, y) + T i{qgq)B (x, y) . (11) 

One may argue that this is not the only way to obtain gauge invariant correlations. In 
particular, the combination (y — x)T^ qgq ^ B (x,y) is the light-cone expression for a gauge 
invariant distribution involving the gluon field strength rather than the covariant derivative. 
On the other hand, these contributions necessarily involve coefficients C{Z gq \ a (x,y) that 
vanish when x = y. Eq.(^) then implies that one can consider instead the combination of 
Eq.(|TTD free of charge. Hence these are the only distributions relevant to our process. Using 
hermitian conjugation, it is easy to show that K\ B is symmetric in x and y and K 2B is 
antisymmetric. We will use these symmetries to simplify our presentation of the coefficient 
functions. 

Let us consider the use of the above approach to the tree diagrams shown in Fig. 3, where 
we have shown only. Those diagrams for Mfj? can be obtained simply by replacing 
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(a) 



(c) 



FIG. 3: Tree contribution to M^ q . The contribution to is obtained by replacing the gluon 
vertex by ij". 



qq 



the gluon interaction it a ^ff_ by the transverse vector vertex Figure 3 corresponds to the 
following on-shell process: a quark and gluon with four- momenta xp^ and (y — x)p' 1 , respec- 
tively, scatter with a photon of four-momentum q^, producing a quark of four-momentum 
yp^ and a forward photon. The Compton amplitude can easily be calculated using the usual 
Feynman rules. When the quark and gluon combine before or after interacting with the 
photons, the intermediate propagator, i/(y drawn with a bar in Fig. 3, is singular and 
needs regularization. Adding an infinitesimal A fl, one arrives at 



i{y p'+X 



i p 1 i fl 
(yp + An) 2 ~ 2A + ~2y 



(12) 



The first term yields zero when acting on the external quark wave function, so only the 
second term contributes. Combining the result from and taking into account cross- 
ing symmetry, we find the following gauge-invariant expression for the tree-level Compton 
amplitude 



x(xb — x) 



[K u (x, y) + K 2i (x, y)\ - (x B -> -x B ) 



(13) 



where if = e 2 — e 2 is the non-singlet part of the quark charge squared and e 2 = J2i e i/ n f- 
The integration has support only for {x, y} G [—1, 1]. When \x B \ < 1, the above expression 
develops an imaginary part through xb — > %b — it- Using the optical theorem, one obtains 
the following tree result for g^ s , 



9t 



NS(0) , 



1 2 r 1 

xb, Q 2 ) = -E e 2 — / dy {K u (x B , y) + K 2i {x B , y)) + (x B -> -x B ) ■ 
* Xb J-i 



^From QCD equations of motion, one can show 

2 f 1 

Aq T (x) = - dy (Kx(x, y) + K 2 (x, y)) , 
x J-i 



where Aqx(x) is defined as 

Aq T (x) 



(14) 



(15) 



(16) 
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This object seems to have a simple physical interpretation. 

The tree level result suggests the following all-order factorization formula for the non- 
singlet part of gr, 

+C2 ^^7, a s ^j i^2i(a:, y)^ + (x B -> -xb) . (17) 

Ci 5 2 can be written in terms of a perturbation series in the strong coupling a s , 

y(») 



The tree result yields 



^ 2 = E^gJ • (is) 



y6(x-x B ) . (19) 




Now we come to the main subject of the paper: one-loop radiative corrections for 
9t S { x b,Q 2 )- All the Feynman diagrams for M^ q are shown in Fig. 4. The first 12 di- 
agrams have color factor Cf, the next 8 come with the factor Cf — Ca/2, and the final 4 
with Ca = N c . Since the external transverse momentum of the quark only goes through 
the quark propagators, the Taylor-expansion of yields the vector vertex insertion on 
the quark line only. These diagrams correspond to the first 20 diagrams in Fig. 4 with the 
external gluon vertex replaced by the vector vertex. An inspection of the one-loop Feynman 
integrals reveals a simple rule: Mj^ is just the Cp part of Mjfg and, hence, the former need 
not be calculated separately. 

Since all Feynman diagrams are computed in bare perturbation theory, the result depends 
on the bare coupling g B . We replace it with a renormalized coupling in the MS scheme, with 
a difference of higher order in a s (n 2 ). After a lengthy calculation involving one-loop integrals 
with up to five internal Feynman propagators, we find 



S? NS (x B ,Q 2 ) = £e, 2 / 



dxdy 



xy 

(1) / x ,) 



M « (JL f JL) KliB ( Xi y) 

\x B x B : 



— )K 2iB (x,y) 
\x B x b 



(x B -> -x B ) . (20) 



Explicit expressions for the M's in the region \x B \ > 1 can be found in the Appendix. 

Now we are ready to show that (x B ,Q 2 ) is factorizable at the one-loop level, i.e., 
the infrared poles 1/e in match the ultraviolet poles in K iB . To this end, we use the 
infrared poles in M; to generate a scale evolution equation for the parton distributions. From 
the result in the Appendix, we find 

d 



dln^L ^ {w^) iKl ^ y ^ 2) + K2iX,y ^ 2)) 
a ^p- j^dxdylcrK^y^ 2 ) 



1 2 2 
+ s 



x(x B — x) x 2 xy 
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2xb 1 xb 
x z y 
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2x, 



y — x \x(y — x) 
+C F K 2 (x,y,iJ 2 
xb — x 



+ 



x(xb 

1 N 



xb - y , 
-1 



X) 



log 



log 



(xb — x 



V x B 
x B -y + x 
xb 



x(xb 



x log 



xb 



x) 
2 



4-2 (l-5? + 



x z 



x z 



X" 



(y-x)(x B -y) V ^b 



4x/ 



2x 



B 



2^/(2/ - x; 
1 N 



y — x \x(y — x) 



+ 



2 



log 
log 



2 (y-x){x B -y) 



( x B - g 
V £b 
xb - y + x 
xb 

^b — y + x 
xb 



log 



(21) 



Expanding the above in the large xb limit, we get the evolution equation for the moments of 
the parton distributions. A detailed check shows that the equation is identical to that in Ref. 
fllCl obtained by studying the ultraviolet divergences present in the twist-three operators. 
In particular, in the large N c limit, the above equation becomes autonomous [|17], |18|| . 

The final step of the calculation is to take the imaginary part of the factorized S^ s (xb, Q 2 ) 
to get a factorized expression for the structure function g^ s (xB,Q 2 ) in the physical region 
xb < 1. With the definition in Eq. (19) we find the following result for the coefficient 
functions, 



£,(1) / x B x b 



C F y { -25{x -x B ) + e[ 1 

,Xb 



3/2 



2x 



B 



y(y 



X) 

2 



+6 



+ 



(xb - x 

y-x 



1/2 

X 

log 



1 

y 



l 



y-x 



(xb - x) 
2xb 1 

X 



_ xb 6xb 

X 2 



xy 



xb 
x 2 



log 



x — x B 
x b 



x b 
2x B 



y-x 



2x 



X 



\y-x 

+5(x - xb) log 



xb - y 
1 - 



j log 

xb 



4x B 

xb - y y-x 
y-x-x B " 



2x 



x B 



x 



x B 
log 



xb 
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xb 



+6(y - x B ) 
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x B 



-e 



y-x 

x B 
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- 1 
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x B 
1 - 
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( 2x 
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xb 
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1 



^log 
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- y (y-x) 



xy 
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x b 



xb 



x 



x B 
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x l 1 -2 lo H— — 

2x d_ii (JL 

(x B — x) V 4 Vxb 



S(y-x B ) 



log 



a; — xb 
x B 



log 



X \ 
Xb) 



(22) 



£,(1) / Xb_ Xb_ 

2 V x ' 2/ , 



C F y <j -25(x - x B ) + ( 1 

v Xb 



x 



3/2 



(x s - x) 



1/2 



x 



1 
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+0 
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x B 

+5(x - x B ) log 
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y-x 
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log 
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Ca 
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log ( — 

xb 
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2x 
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log 
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log 
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log 



J/ 

x — xb 
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y-x 



- 1 
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log 



X 

xb 



log 



X 

xb 



(23) 



The definition of the + functions can be found in Ref. ||. The support for the parton 
correlations limits y — x to the interval [—1, 1]. 

To check the Burkhardt-Cottingham sum rule |0|, we integrate g^ (xb,Q 2 ) over xb- 
Assuming the integration over x and y can be interchanged with that of xb, one obtains 

f Q dx B g» s {x B ,Q 2 ) ^E^ 1 - \ Cf ^) (PS$n±-wl>i\PS) ■ (24) 

Here the coefficient 7/2 reduces to 3/2 if we define 75 so that the non-singlet axial current is 
conserved. Compared with the factorization formula for gi(x B , Q 2 ), we have the Burkhardt- 
Cottingham sum rule at one loop 



f 1 dx B g? s (x B , Q 2 ) = . 
J 



(25) 



If the order of integration cannot be interchanged because of the singular behevior of the 
parton distributions at small x and y, the above sum rule may be violated. Indeed some 
small xb study does indicate such singular behavior [2~0|. 

Finally, we consider the next-to-leading order correction to the non-singlet part of the x 2 
moment of gr(x B , Q 2 )- in the leading order, it is well known : 



jf 1 dx xVHx, Q 2 ) = ~ £ e 2 Qa 2l (Q 2 ) + d 2l {Q' 



(26) 
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where a 2 is the second moment of the g\(x,Q 2 ) structure function and d 2 is a twist-three 
matrix element [13]. Using the coefficient functions found above, we see that 



I 1 dxx 2 g" s (x,Q 7 
Jo 




1 + 



a s (Q 2 ) 7 



21 



1 + 



4tt 12 

a s (Q 2 



C h 



4lT 



27„ 29 

T c A - -C F 



(27) 



Using the next-to-leading result for gi(x, Q 2 )0, we find 
jjxx? (j? s (*,Q 2 ) - i 9 f s (x,Q 2 )) = jE^a (l + 



-g!l(f C .-f Cp )]•(-) 



Notice that the combination of #t and g\ relevant to d 2 receives no radiative correction. 

To summarize, we have extended the leading-twist factorization formalism to higher twist. 
The extension mainly involves a correct identification of the intermediate parton states and 
arranging different calculations into gauge invariant combinations. Because of the multipar- 
tons in the initial state, higher-twist perturbative calculations are, in general, much more 
complicated than the leading twist cases. As an example, we have presented the one-loop 
coefficient functions the twist-three process involving forward Compton scattering on a trans- 
versely polarized nucleon target. This result can be combined with a two- loop calculation 
of the twist-three evolution equation to yield a next-to-leading order formalism for the g 2 
structure function. 
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APPENDIX 

In this appendix, we present the tree and one-loop Compton amplitudes for virtual photon 
scattering on a transversely polarized nucleon in the unfactorized form. To simplify the 
expressions, we assume \xb\ > 1 so that the amplitudes are purely real. With e 0123 = +1, 
we write the Compton tensor as 



T^ = -i€^ aP q a S ± p-S T (xB,Ql 2 



(29) 



In QCD, we can write, 

S? s (x B , Q 2 ) 



t J -1 xy L \x B XbJ 

+M 2 (— ,-9- )K 2iB (x,y) 
\x B x B J 



- (xb -> -x B ) 



(30) 
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where the ICs are unrenormalized parton distributions as defined in the text, M's are 
perturbation series in a s and have infrared poles. 
At tree level, the result is well known, 



Mi \x,y) = M?>(x,y) = 



(0), 



Jy_ 

l-x 



(31) 



At the one-loop level, the following amplitudes are associated with distributions with a 
covariant derivative, 



x log(l - x) + 
-2 



?Y /2 [2 



-y 



(l-x) 
( 2 



x 
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[ e 
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+ — + 
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log(l — x) 
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H I log(l - x) 
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1 - 



2x 
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log(l -y + x) 
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(32) 
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+— ^ 1 + log 1 - y + x) 

y-x \ l-yj 
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(i/ - ^X 1 - v) 

where e = 4 — d, p, 2 = 47re _TB yU 2 , and 7^ is the Euler constant. 
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